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NMR Evidence for Gapped Spin Excitations in Metallic Carbon Nanotubes 
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We report on the spin dynamics of 13 C isotope enriched inner-walls in double-wall carbon nan- 
otubes (DWCNT) using 13 C nuclear magnetic resonance (NMR). Contrary to expectations, we find 
that our data set implies that the spin- lattice relaxation time (Ti) has the same temperature (T) 
and magnetic field (H) dependence for most of the innerwall nanotubes detected by NMR. In the 
high temperature regime (T > 150 K), we find that the T and H dependence of 1/TiT is consistent 
with a ID metallic chain. For T < 150 K we find a significant increase in 1/TiT with decreasing T, 
followed by a sharp drop below ~ 20 K. The data clearly indicates the formation of a gap in the 
spin excitation spectrum, where the gap value 2A ~ 40 K (= 3.7 meV) is H independent. 

PACS numbers: 71.20.Tx, 61.46.+W, 73.22.-f, 76.60.-k 
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The electronic properties of carbon nanotubes have 
been of a topic of intense investigation ever since their 
discovery in early 1990's. According to band-structure 
calculations the basic electronic structure of single-wall 
carbon nanotubes (SWCNT) is expected to depend on 
the chiral wrapping vector (n, m) across the graphene 
plane, where tubes for which (2n + m)/3 = integer 
are metallic, while all other tubes are semiconducting 
00 HQ with a large ~1 eV gap @. While STM and 
transport measurements on isolated tubes demonstrate 
the diversity of tube properties, significant measurements 
on macroscopic amounts of tubes are only possible in se- 
lected cases. Photoemission spectroscopy (PES) mea- 
surements on metallic tubes in bundles 0, Q suggest 
that strong electron-electron correlations can lead to a 
Tomonaga-Luttinger-liquid (TLL) state. Recently, dou- 
ble wall carbon nanotubes (DWCNT) have been syn- 
thesized by filling SWCNT with fullerenes (so called 
"peapods" [8j) followed by a high temperature reaction 
which merges the fullerenes into inner tubes 0, llfj| . 
These DWCNT have some exceptional properties since 
the inner tubes are accommodated in a highly shielded 
environment under clean room conditions |ll| . Raman 
experiments performed even on bucky paper material al- 
lows one to detect some significant properties of these 
inner tubes due to their small diameter (high curvature) . 

Nuclear magnetic resonance (NMR) is usually an ex- 
cellent technique for probing the electronic properties at 
the Fermi level of metallic systems, take for instance 
conducting polymers, fullerenes, and high temperature 
superconductors to mention a few. However the 1.1% 
natural abundance of 13 C with nuclear spin 1=1/2 lim- 
its the sensitivity of such experiments. Data taken on 
SWCNT essentially evidence a large distribution of prop- 
erties since samples of identical tubes are presently out 
of reach. In this report, selective enrichment of the in- 
ner shells using 13 C isotope enriched fullerenes 0, El 
in the "peapod" synthesis route is used to probe the spe- 
cific properties of the inner tubes. The 13 C enrichment 



allows us to increase the 13 C NMR sensitivity by two or- 
ders of magnitude, and furthermore achieve selective 13 C 
enrichment of the inner shells alone. This provides us 
with the possibility of singling out the electronic prop- 
erties of these inner tubes for the first time. We show 
that, although these tubes are distributed in diameter 
and chirality, their electronic properties display a strik- 
ingly homogeneous behaviour. The magnetic properties 
of these inner wall nanotubes behave as for a ID metal 
at room T, but exhibits a pronounced gap below ~20 K. 
This unexpected result reveals that this specific macro- 
scopic collection of carbon nanotubes is an object display- 
ing original physical properties worth studying in more 
detail with macroscopic experimental techniques. 

All 13 C NMR data in this report were taken with the 
sample sealed in a 6 mm diameter glass tube filled with 
200 mbar of high purity Helium gas. Details of the syn- 
thesis techniques and various independent experimental 
evidences for the formation of 89% 13 C isotope enriched 
inner-walls inside natural 1.1% 13 C enriched outer walls 
are reported elsewhere . We probed the low frequency 
spin dynamics (or low energy spin excitations, equiva- 
lently) of the inner-tubes using the spin lattice relaxation 
time, T\, defined as the characteristic time it takes the 
13 C nuclear magnetization to recover after saturation. 
The signal intensity after saturation, S(t), was deduced 
by integrating the fast Fourier transform of half the spin- 
echo for different delay times t. All data were taken with 
excitation pulse lengths it/2 = 3.0/xs and short pulse sep- 
aration times of r = 15/xs |l4|. We obtained the value 
of Ti by fitting the t dependence of S(t) to the form 
S{t) = S a - S b - M(t), where S a ^ S b (> 0) are arbi- 
trary signal amplitudes, and 



M(t) = cxp - (i/Tf) 



(1) 



is the reduced magnetization recovery of the 13 C nuclear 
spins. Fig. ^ shows the results of M(t) for the inner- 
tubes as a function of the scaled delay time t/Tf, un- 
der various experimental conditions listed in the Figure. 
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We find that M(t) does not follow the single exponential 
form with (3—1 (dashed line), but instead fits well to 
the stretched exponential form with f3 ~ 0.65(5) which 
implies a distribution in underlying relaxation times Ti 
across the sample. In such cases, Tf in Eq. is directly 
proportional to the mean value T\ of the T± distribution 
as such Tf = Tj • f3/T(l/(3), where T is the gamma func- 
tion. We display the data in Fig. Qon a semi-log scale for 
the time axis in order to accentuate the data for earlier 
decay times and to illustrate the collapse of the data set 
for the upper 90% of the NMR signal. We find that the 
upper 90% of the M(t) data is consistent with constant 
ft ~ 0.65(5) (see inset), implying a constant underlying 
distribution in T\ for a large range of experimental condi- 
tions. The lower 10% of the M(t) data (corresponding to 
longer delay times) comes from the non-enriched outer- 
walls which, as a result of their larger diameters, have 
much longer relaxation times under similar experimental 
conditions [3 03, El E3 • 

Two distinct origins for the multi-exponential magneti- 
zation recovery can be considered. The first is due to the 
powder average over the spatial anisotropy in T\. The 
distribution is independent of the tube properties, and 
can also be found in the 13 C NMR data for alkali doped 
fullerenes A„C6o U3Ei3| Given the similar diameter of 
C-60 (d = 0.71 nm) to the average inner-wall diameter 
(d = 0.7 nm |12i in this report, we can expect com- 
parable bonding effects for the electron orbitals. It has 
been shown that in A„C6o the Xi for 13 C is dominated 
by dipole-dipole interactions between the electron spin 
in the ppir bond and the 13 C nuclear spin [2lj| . In this 
case, the relaxation depends on the orientation of the p^ 
orbital (which is perpendicular to the tube surface) and 
the external magnetic field, and therefore contributes to 
the multi-exponential form of magnetization recovery for 
a powder average. This resultant Ti distribution is inde- 
pendent of T and H. 

Another source of multi-exponential recovery is from a 
distribution of the inner tube properties themselves, such 
as their diameter. According to Raman scattering, the 
inner tubes have a mean diameter of d ~ 0.7 nm with a 
standard deviation of a ~ 0.1 nm Within this 

distribution lies a variety of tubes with different chiral- 
ity and one can a priori expect to find metallic as well 
as semiconducting tubes . If both semiconducting and 
metallic inner-tubes existed in our sample, one would ex- 
pect the ratio of the Ti's between the different tubes to 
increase exponentially with decreasing T below the semi- 
conducting gap (~ 5000 K 0), which would drastically 
change the underlying Ti distribution with decreasing 
T. This change would manifest itself as a large change 
in the shape of the recovery M(t), however, as shown 
in Fig. n this is not the case. We can therefore rule 
out the possibility of two components in T\ with differ- 
ent T dependences, and instead we conclude that all Ti 
components exhibit the same T and H dependence within 



experimental scattering. 

The Ti distribution in the sample, whether it arises 
from anisotropy or diameter variations (or both), shows 
a uniform T and H dependence. It is therefore appro- 
priate to follow the T and H dependence of the mean 
value of the distribution (T-f in Eq. JIJ), and thereby 
get insight into the homogenous electronic state of the 
inner tubes. In order to avoid unnecessary experimental 
scattering in Tf, we then go back and fit all the M(t) 
data to Eq. with a fixed value of j3 = 0.65. We plot 
the resulting temperature dependence of 1/TfT in Fig. 
[5] for two different values of the magnetic field H . We 
can immediately separate the data into two temperature 
regimes; the high temperature regime > 150 K, and the 
low T regime < 150 K. At high temperatures we find 
that 1/ TfT is independent of T which indicates a metal- 
lic state |14j. which given the arguments above implies 
that all of the inner tubes are metallic. We also find a 
strong field dependence for T\ which is best illustrated 
by plotting the high temperature value of 1 /Tf T against 
l/\/~H for H values ranging from 1.2 Tesla to 9.3 Tesla, 
as shown in Fig. |3| We find that the data fit well to the 
form 

_L = A + B-$=, (2) 

TfT ViT 1 ; 

where A and B are constants, which is very sug gestive of 
a ID spin diffusion mechanism for Ti jHHlliJ. B/\/~H 
corresponds to the diffusive contribution to the relaxation 
originating from the long wavelength (i.e. q ~ 0) modes, 
while A corresponds to the non-diffusive contributions 
from q > modes. A cutoff to the divergence in Eq. @ 
as H — > is often encountered in ID spin chain systems 
[22I I23I IzS ] due to inter-chain coupling. In the present 
case, we can postulate that electron tunneling between 
inner to outer walls could cause similar cutoff effects, 
however the data down to the lowest field of H=1.2 Tesla 
indicates that the cutoff has not been reached yet. We 
therefore conclude that the high-temperature regime is 
consistent with a ID metallic chain. 

The origin of the unusual T dependence of 1/TfT in 
the low temperature regime (< 150 K) is not immediately 
obvious. We can however rule out certain possibilities. 
Firstly, we can rule out the possibility of an activation 
type mechanism where T± is dominated by fluctuating 
hyperfine fields which are slowing down with decreasing 
T. If this were the case the temperature where 1/TfT 
reached its maximum would shift with the resonance fre- 
quency u> E3> or the applied magnetic field uj — j n H 
(7 n /27r = 10.71 MHz/Tesla is the gyromagnetic ratio for 
13 C), equivalently. As shown in Fig. |2 however, we 
find no evidence of a shift in the peak temperature with 
H. Furthermore, at low temperatures 1/TfT is found to 
drop below its high temperature value which rules out 
the possibility of an activation contribution plus a T in- 
dependent contribution. Secondly, we also rule out the 
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possibility that the T dependence of 1 /TfT is a result of 
paramagnetic centers which can arise from wall defects or 
impurity spins. The fact that a pronounced gap exists in 
1/TfT implies a pronounced gap in the low energy spin 
excitation spectrum, which cannot be explained by the 
presence of paramagnetic centers. We note that at the 
lowest temperatures < 5 K (not shown), Tf becomes so 
long (> 300 s) that the low energy spin excitations spe- 
cific to the homogeneous properties of the inner-tubes 
become inefficient, and other excitations take over, pos- 
sibly defect related. In such cases we find that the shape 
of M{t) is no longer universal and that the underlying 
distribution in T\ is smoothed out, possibly as a result of 
nuclear spin-diffusion. 

Having ruled out the above possibilities, we are then 
lead to consider the simplest explanation for the exper- 
imental data using a non-interacting electron model of 
a ID semiconductor with a small secondary gap (SG). 
The SG may be a result of the finite inner- wall curva- 
ture Lji la 12.4 l'2("il| , or perhaps the applied magnetic field 
itself [2jJ . We can fit the 1/TfT data using this non- 
interacting model with only one free parameter, the ho- 
mogeneous SG, 2A. We start by taking the normalized 
form of the gapped ID density-of-states n{E) 

n(E) = fOT| ^ l>A (3) 

I otherwise 

also known as the van-Hove singularity (E is taken with 
respect to the Fermi energy). We then use Eq. J2J to 
calculate 1/TfT as such 

= a(w) J" n(E)n(E + u) dE, (4) 

where E and u> are in temperature units for clarity, / 
is the Fermi function / = [exp(£'/T) + 1] , and the 
amplitude factor a(tu) — A + B' j sfw is taken directly 
from Eq. J3J> and Fig. E](where B' = 4.53- lO^s" 1 !^ 1 / 2 
for uj = H in temperature units). We note that factoring 
out the diffusion effects from the integral in Eq. is an 
approximation valid only if A and B are T independent. 
Eq. 10} cannot be solved analytically, therefore we resort 
to numerical integration. The results of the best fit to Eq. 
Q are presented in Fig. [3 where 2A = 43(3) K (= 3.7 
meV) is found to be H independent (within experimental 
scattering) between 9.3 and 3.6 Tesla. We note that at 
the largest external field of 9.3 Tesla, ui = 4.5 mK <C 
A,T, however, oj must be retained inside the integral. 
This is a consequence of the one dimensionality which 
yields a logarithmic divergence inside the integral of the 
form ln(T/w) for T < A. 

What could possibly be the origin of the observed gap! 
Tight binding calculations predict that applied magnetic 
fields can induce SG's of similar magnitude for metallic 
SWCNT [23. However, such a scenario is excluded here 
from the absence of field dependence of the observed gap. 



Our data would be more consistent with a curvature in- 
duced SG for metallic tubes 0, 0, |2^, |2(|, however for 
our typical inner-tubes the predicted values, ~ 100 meV, 
are over an order of magnitude larger than our exper- 
imental data. Other scenarios, such as quantization of 
levels due to finite short lengths of the nanotubes could 
be considered as well, however, in all these cases a be- 
haviour independent of tube size and chirality is certainly 
not expected. 

This leads us to consider the effect of electron-electron 
interactions for the metallic inner tubes. It has been 
predicted that electron-electron correlations and a TLL 
state leads to an increase in \jT\T with decreasing T [2^ , 
which is a direct consequence of the ID electronic state. 
The correlated ID nature may also lead to a Peierls in- 
stability Q with the opening of a small collective gap 2A 
and a sharp drop in \ jT\T below A ~ 20 K. Therefore, 
the presence of both a TLL state and a Peierls instability 
could possibly account for the data, although here again, 
the independence on tube geometry should be accounted 
for. 

In conclusion, we have shown that the T\ recovery 
data indicate that most of the inner-tubes have similar 
T and H dependences, with no indication of a metal- 
lic/semiconductor separation due to chirality distribu- 
tions. At high temperatures (T > 150 K) 1/TfT of the 
inner tubes exhibit a metallic ID spin diffusion state, 
with no low- field cutoff down to 1.2 Tesla. This metal- 
licity could result from charge transfer from the outer 
to the inner tubes, however this speculation ought to be 
confirmed by independent experiments and theoretical 
calculations. Below ^150 K, 1/TfT increases dramati- 
cally with decreasing T, and a gap in the spin excitation 
spectrum is found below A ~ 20 K. We list various inter- 
pretations for this temperature dependence, ranging from 
a non-interacting secondary band-gap model to a ID cor- 
related electron model with a collective gap (possibly a 
Peierls instability). Firstly, these results should stimu- 
late further experimental investigations on diversely syn- 
thesized DWCNT in order to check whether these ob- 
servations are specific to the "peapod" synthesis route. 
Secondly, theoretical work on the incidence of ID cor- 
relation effects for inner-wall nanotubes inside DWCNT 
should be helpful in sorting out the origin of our aston- 
ishing experimental evidence. 
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FIG. 1: Reduced nuclear magnetization recovery, M(t), as a 
function of the scaled delay time i/Tf (see Eq. for various 
experimental conditions. Both axes are dimensionless. Solid 
grey curve shows stretched exponential fit with (3 — 0.65, 
while grey dashed curve shows single exponential with (3 = 1. 
Inset shows temperature dependence of the best fit values of 
(3 at 3.6 Tesla (•) and 9.3 Tesla (o), and average value of the 
data set (3 = 0.65 (solid line). 



FIG. 2: Temperature dependence of spin-lattice relaxation 
rate divided by temperature, 1/TfT, in units of (10 3 x s _1 K 
Grey curves are best fits to Eq. 0} with 2A = 46.8(40.2) 
K for H = 3.6(9.3) Tesla, respectively. 

FIG. 3: 1/TfT, in units of (10 3 x s" 1 K ~ 1 ), at fixed T = 290 
K, plotted as a function of 1/s/H, in units of (Tesla -1 / 2 ). 
Linear fit corresponds to 1/TfT = A + B/VH with B = 
0.00206 (Tesla 1/2 s" 1 KT 1 ) and A = 0.00028 (s -1 KT 1 ). 
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